Plasmons in the presence of Tamm-Shockley states with Rashba sphtting at noble 

metal surfaces 
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Au(lll) or similar noble metal surfaces feature Tamm-Shockley surface states that are known to 
possess considerable spin-orbit splitting of the Rashba type of order A = 0.1 eV. When interacting 
with an electromagnetic field such states are expected to have resonances when the frequency of 
the field is near the energy of the spin-orbit splitting A. They originate from the intersubband 
transitions between spin-split subbands and can be observed in the frequency dependence of the 
surface impedance. Plasmons in thin metal films are gapless and can be strongly affected by these 
spin resonances, acquiring significant modification of the spectrum when it intersects the = A 
line. Finally, an interesting demonstration of the intersubband resonances can be achieved when 
metal films are coated with ionic dielectrics that have a frequency of longitudinal/transverse optical 
phonons above/below A. The dielectric function between the two optical phonon frequencies is 
negative which forbids propagation of conventional plasmon-polaritons. However, the presence of 
spin-orbit-split surface states allows plasmon-polaritons to exist in this otherwise forbidden range 
of frequencies. 
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PACS numbers: 72.25.-b, 73.20.-r, 73.50.Mx, 78.66.B2 



I. INTRODUCTION 

Nanoplasmonics is a novel field that emerged at the 
confluence of optics and condensed matter physicsii^i^. 
Its ultimate goal is the development of high resolu- 
tion imaging methods by means of plasmon-enhanced 
near-field optical measurements. Plasmons are collective 
charge excitations of electron liquids that are induced by 
external electric fields. In a hulk metal they are purely 
longitudinal and have a gapped spectrum given by the 
well-known Langmuir expression^, 



(1) 



in terms of bulk electron density A'' and mass m. At 
the interface between a metal and an insulator (or vac- 
uum) a distinct but closely related excitation can prop- 
agate - a surface plasmon, whose frequency is reduced 
significantly, uig — fl/ + k, by the dielectric constant 
of the insulator i^^. In thin metallic films such surface 
plasmons propagating near opposite interfaces become 
hybridized and split into a symmetric mode (with oscil- 
lating charges at the interfaces having the same sign) and 
anti-symmetric mode (opposite sign). As a result the fre- 
quency of the anti-symmetric mode increases while the 
symmetric mode becomes gapless. 

An electron liquid confined to two dimensions (2DEG), 
created in semiconductor heterostructures, features a 
plasmon spectrum remarkably different from its three- 
dimensional analog. In particular, due to weaker screen- 
ing the plasmon spectrum of a single 2D layer is gapless^, 



k. 



(2) 



E(eV) 




0.8 1 1.2 
p (Angstrom ^) 

FIG. 1: (color online) Energy spectrum of Au(lll) films rel- 
ative to the Fermi level (E=0): solid lines show the spectra of 
spin-orbit-split Tamm-Shockley states, E = /2m±ap/2 vs. 
the in-plane momentum. Dashed lines show the bulk electron 
spectrum E = (/ +pl)/2m for p^ = 0;0.2A"\ 



and TO* is its effective mass, typically significantly lower 
than the vacuum electron mass toq. Another celebrated 
property of a two-dimensional electron system is the spin- 
orbit interaction. As known from the relativistic Dirac 
equation, in the presence of any potential U{r) the elec- 
tron spin operator s is coupled to its momentum p via 
the Hamiltonian, 



Hso = -AVJ7 • (s X p). 



(3) 



where k is a dielectric constant of a host semiconduc- 
tor, n is the planar density of the 2D electron system 



Hereinafter we use units with h = 1. 

In a bulk crystal an aperiodic part of the potential 
U originates from impurities and is relatively weak in a 
clean crystal. In a 2DEG, on the other hand, it is in- 
trinsically present as it comes from a confining interface 
potential which creates the 2DEG*^. Interestingly, the 
value of the coupling constant A in a typical GaAs semi- 
conductor is six orders of magnitude stronger than in the 
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vacuum (and of opposite sign), where it is found from the 
Dirac equation to be A = — 1/2tooC^ ~ — 2x10^^ A^. De- 
spite this favorable fact, the spin-orbit splitting of a 2D 
electron spectrum is rather weak^: at the Fermi level, 
A ~ 0.1 — 1 meV, since the typical Fermi momenta are 
small, ~ 0.01 A^^. As a result, the spin-orbit splitting 
leads to relatively weak corrections^"'^^ to the 2D plas- 
mon spectrum, Eq. and is difficult to observe exper- 
imentally. 

In the present paper we report on a novel phenomenon, 
which combines the properties of bulk metals and two- 
dimensional electron systems, and predict a resonant cou- 
pling of individual electron spin degrees of freedom to 
a collective charge excitation originating from the spin- 
orbit interaction. 

Photoemission experiments performed on (111) sur- 
faces of goldi^^ii^'^ii^ reveal the existence of L-gap 
surface Tamm-Shockleyiii^ states that constitute a two- 
dimensional electron gas. Its properties are described by 
the Fermi momentum pp « 0.2 A'^^ and effective mass 
m* = 0.2mcr^; for comparison, the Fermi momentum for 
bulk electrons is w 1.2 A~^. Of significant importance to 
us is the observation that the spin-orbit splitting of the 
surface states is unusually strong, which is measured to 
be A « 0.1 eV at the Fermi level. Recent experiments^^ 
show that Bi/Ag(lll) surface alloys can have a splitting 
which reaches A w 0.2 eV. As we demonstrate below it is 
this strong spin-orbit interaction that leads to coupling of 
electron spin to charge oscillations via resonant electron 
transitions between spin-split subbands. Fig. 1 shows 
the energy spectrum of both bulk and surface states for 
Au(lll). 

A. Geometry of the system 

We will consider two geometries: a) a noble metal 
film of thickness d embedded into a insulator with a di- 
electric function k{lu); b) a single metal interface with 
a vacuum. For the former geometry we are interested 
in the plasmon-polariton spectrum, while the geome- 
try b) is discussed in relation to the observation of the 
surface state resonances in the reflection coefficient (or 
impedance) of electromagnetic radiation incident on the 
surface. 

Both geometries require first that we analyze the dy- 
namic response of both surface and bulk electrons. The 
latter are conventionally described by the Drude dielec- 
tric function 



where the bulk plasma frequency in gold is = 1.3 x 
lO^^s"^. The electron momentum relaxation rate l/r 
is typically much smaller and depends on temperature 
(phonon scattering) and crystal quality (impurities) . Ex- 
pression Eq. ^ is the approximation that neglects vac- 
uum contribution and d-band electrons. However, the 



resonant phenomena addressed in our paper occur at 
infrared frequencies, to ^ 100 meV. At these frequen- 
cies the s-band contribution is large and dominant, and 
Eq. (dl to be a good approximation. 
The Hamiltonian of the surface states 

-H^ = ^— ^ + a(s X p)^, p^{p^,py), (5) 

describes the in-plane dynamics of confined electrons, 
with the effective spin-orbit coupling constant given by 
the average of Eq. ([3]) over the direction perpendicular 
to the surface {z), 

oo 

a^-Ay dz\Mz)\^^, (6) 

— oo 

where ijjp is the wave function of a surface state. We 
make an approximation (consistent with the photoemis- 
sion data) that a is momentum-independent. The nu- 
merical value of a is related to the observed value of 
spin-orbit splitting according to A = app- 



II. DYNAMIC RESPONSE OF SURFACE 
STATES 

Qualitatively, the interaction of the electromagnetic 
field (plasmon, infrared beam reflected from the surface, 
etc.) with our system can be understood as follows. An 
electric field along the x-direction causes charge oscilla- 
tions of both bulk and surface electrons and an induced 
electric current jx- It is the resonant contribution to this 
current from the interband 2D transitions which is of the 
most interest to us here. The result of the appearance of 
the electric current is two-fold: 

i) According to Ampere's law, a magnetic field along 
the y-direction appears that leads to the interface Zee- 
man magnetization My. 

ii) As seen from Eq. ([5]) a drift along the x-direction can 
be considered as effective "Zeeman" field pointed along 
the y-direction. Consequentially, it induces additional 
net electron spin polarization and, hence, magnetization 

Subsequently, the surface magnetization and the ac- 
companying ac magnetization induce an additional ac 
electric field (along the x-direction) via Faraday's law 
which acts upon both bulk and surface electrons. It is 
crucial that all these effects are resonantly enhanced for 
frequencies near the frequency of inter-subband transi- 
tions, a; ~ A. 

Let us now address the problem of dynamic response 
of 2D electron states quantitatively. Most simply it can 
be done by means of a kinetic equation. Due to the spin 
structure of the Hamiltonian ([5]), the electron distribu- 
tion function /p is a 2 x 2 matrix in spin space. The 
corresponding equation for /p in the presence of both ac 
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electric and magnetic fields has the form Ref . |2l|, 



dfp 

dt 



[finjp] = ^gfJ'B[^yJp]Hy 



ei?.^, (7) 

OPx 



where Ap = ap, fjn — ny&x — rixCTy, with n = p/p being 
the direction of electron momentum and it = 2s the set 
of Pauli matrices; /is = \e\/2mc is the Bohr magneton 
and g is the gyromagnetic ratio of surface electrons. 

A formal solution to the equation ([7]) in the frequency 
domain reads (see Appendix A for more details) 



(8) 



where ICp denotes the right-hand side of Eq. ^ . To the 
linear order in the fields Hy and E.^ it is sufficient to 
utilize the equilibrium distribution function in the right- 
hand side of Eq. /p"^ = ^ '>ln)n±, with n± — 
"■f(p^/2to*± Ap — /i) denoting the Fcrmi-Dirac functions 
for the two spin-split subbands. 

Using Eq. ([8]) it is now starightforward to calculate the 

surface magnetization M = — ^.g/isTr J ^^J^2 fpO', and 

current density jx = e Tr/ j§:;^fp{^ - ^aay). Note 
that the electron velocity contains spin operator. After 
momentum integration we obtain for the 2D electric cur- 
rent and magnetization, 



jx = <j{cu)Ex{0) - iojP{oj)Hy{0), 
My^P{u;)Ex{0) + xi^)Hy{0). 



(9) 



Here we explicitly emphasized that the fields are to be 
taken at the interface {z = 0) (see also the discussion 
in the next section.) Here the electric conductivity of 
surface states (A — Ap^), 



+ 



(10) 



where n = p|n/27r is the density of 2D electrons, consists 
of the usual intrasubband Drude conductivity (first term) 
and the contribution of resonant intersubband transi- 
tions (second term). The dynamic magnetic susceptibil- 
ity x(a-') also features a similar resonant structure, which 
with the help of Eq. ([8]) is found to be 



X(w) = g Ms 



A2 



StT - A2 ■ 



Finally, the cross-susceptibility has the form, 
iem* gfiBCtA^ 



Stt u;(w2_A2)- 



(11) 



(12) 



Note that the off-diagonal terms in Eqs. ([9]) are related 
to each other, in agreement with the Onsager theorem. 

In the presence of two surfaces, expressions similar to 
the above Eqs. ([9lfT2|l hold for the other surface as well 
(with the reversal of the sign of a, which corresponds to 
the reversal of the sign of dU jdz in Eq. ^ . 
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FIG. 2: (color online) Frequency dependence of the abso- 
lute value of the surface Au(lll) impedance \Z(lij) j Z{i{u])\, 
as related to the impedance Zo{uj) = 1/ —e{u]) in the ab- 
sence of Tamm-Shockley states, for different values of the elec- 

(blue dashed line), 



s 



tron scattering rate: 1/r = 1 x 10 

3 X lO^^s"^ (red dotted line), 1 x 10"s"^ (black solid line). 
The energy of spin-orbit-induced intersubband resonance is 
A = lOOmeV = 1.5 x 10"s~\ 



III. IMPEDANCE OF Au(lll) SURFACE 

The most straightforward way to observe resonances, 
Eqs. PTHl - (fT^ . associated with the intersubband tran- 
sitions is to measure the reflection of the electromag- 
netic wave incident from vacuum on a surface of a bulk 
metal. Since typical photon momenta A/c is consider- 
ably smaller than the characteristic 2D electron momen- 
tum scales of the problem, pp and A/vp, it is sufficient 
to consider a case of normal incidence and refiection. As- 
suming the incident wave has frequency uj and amplitude 
Eq, we can write the electric field as 



Ex{z) 



z/c 



EqC " -r £J2 

El exp (y>z^J~^l c) 



< z, 
z < 0. 



(13) 



Accordingly, the magnetic field Ely is found from 
Maxwell's equations. Inside the metal it reads dHy/dz — 
iojsEx/c] and similarly in the vacuum (where e 1). As 
a result. 



Hy{z) = 



-EQe^^^""!" ^ i;2e*"^/^ < z, 
—i\f^^E\ exp {ujz^'-e/c), z < 0. 



(14) 



Boundary conditions relate the discontinuities of the 
fields Ex and Hy to the surface current and magneti- 
zation given by Eqs. ([9]), 




El 



En — E2 — 



(15) 
(16) 



The right-hand side of the boundary conditions psp are 
not well defined as the surface electric current and mag- 
netization, Eqs. ([H]), are given in terms of electric Ex{Q) 
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and magnetic Hy{0) fields right at the surface. The lat- 
ter, however, are discontinuous at z = 0. This means 
that the response of surface electrons has to be in princi- 
ple found from a solution of a three-dimensional problem 
which accounts for the z-dcpendence of the electron den- 
sity distribution, dynamics of surface states, as well as 
Maxwell's equations in a self-consistent way. Such anal- 
ysis is well beyond the scope of the present paper and 
should be a subject of a separate work. In this paper 
we adopt a phenomenological approach and approximate 
that the 2D electron density is constant across the inter- 
face so that that the effective surface electric and mag- 
netic fields in Eqs. ([9]) are given by the corresponding 
mean arithmetic values of the magnitude of the fields 
immediately above and below the surface. 



E40) 



Eq + El + E2 



Hy{0) 



E2-E0- iV^Ei 



In Appendix B we show how these boundary conditions 
follow from the approximation of 2D electron density con- 
stant across the interface. Equations p5)) together with 
([9|) and (fT7|l now give a system of two coupled linear 
equations for Ei and £2- The experimentally measur- 
able quantity, surface impedance Z(lu), is given by the 
ratio of the electric and magnetic fields aX z — 0+, 



Eq + E2 



HyiO+) Eo 
After simple calculation we obtain 



E2 



(18) 



Z{lo) 



l + (^/3)' + ^x(v^-T) 



1 



Substituting now the values of cr, x ^-nd /3 found in the 
preceding section we observe that the doubly-resonant 
terms, oc are exactly canceled by the corresponding 
contributions from the terms oc ax- Of the remaining 
terms only the cr-term in the denominator is to be kept, 
as the x-terms are small by virtue of {vf/c)'^{^/^) ^ 1, 



Z{io) 



1 



(20) 



The resonant feature in Eq. ([^0]) is rather narrow though 
sharp. A finite electron scattering rate leads to its broad- 
ening. 



A. Electron scattering 

The scattering of surface electrons off phonons, impu- 
rities and surface roughness can be be accounted for by 
a collision integral in the kinetic equation ([7]). In a re- 
laxation time approximation this yields the substitution 
u! uj + i/r in the response functions (fTUl) - P^ . The 
broadening of the bulk electron response is accounted for 



by the imaginary part in the dielectric function ([4]). For 
estimates we will assume that the scattering rate l/r is 
the same for bulk and surface states. Fig. [2] illustrates 
the dependence of the absolute value of the impedance 
|Z(ti>)| close to the resonant frequency for different values 
of the scattering rate. Typical scattering rates in noble 
metals at room temperature are ~ 10"'^'' s~^. In order to 
achieve the desired scattering rate of order ^ 10^^ s~^ one 
should make measurements at low, T < 10 K, tempera- 
tures. Due to the screening from bulk electrons, which 
gives large values for the refraction index w i^/A, 
the relative correction to the surface impedance is small 
but within the capabilities of modern optical detection 
methods. 

An alternative setup for the observation of intersub- 
band surface resonances is the measurement of the trans- 
mission coefficient through a thin metal film, rather than 
of a reflection (impedance) from a bulk metal. The corre- 
sponding calculations can be performed in complete anal- 
ogy to the above analysis. 

Finally, we note that the applicability of the linear re- 
sponse close to the resonance requires that the first-order 
variation of the electronic distribution function is small 
compared with the equilibrium Fermi-Dirac distribution. 
Estimating the right-hand side of Eq. ([S]) we obtain that 
close to the resonance the following condition should be 
satisfied, eEvp/^ max(|a; — A|, l/r). This condi- 
tion has a simple physical meaning - the energy acquired 
over the spin-orbital distance vp/A in the electric field 
should be smaller than the detuning to the resonance or 
the scattering rate. 



IV. PLASMON SPECTRUM OF A METAL FILM 

We now consider the role of intersubband resonances 
on the properties of plasmons in thin Au(lll) films, em- 
bedded into a dielectric, see Fig. [31 The plasmon is prop- 
agating along the a;-direction. In general, it is accompa- 
nied by Ex, Ez components of electric field as well as an 
Hy component of magnetic field. As illustrated in Fig. [3] 
there are two modes in a film, corresponding to symmet- 
ric and anti-symmetric alignment of the electric field near 
the two surfaces. 

Let us begin by estimating the relative importance of 
the two resonant contributions, arising from the surface 
electric current and surface magnetization, cf. Eqs. ([9]). 
Integrating Faraday's law, V x E = ia;(H -I- 47rM)/c, 
across an interface, we can obtain that the additional 
electric field induced by the oscillating surface magnetiza- 
tion is JE'i^^ ~ ojMy/c. Similarly, according to Ampere's 
law the additional contribution to the magnetic field due 
to the surface current is 5Hy ~ Sjx/c, where 5jx of inter- 
est is the second (resonant) term in Eq. ([10]). Consequen- 
tially, this magnetic field yields an additional electric field 
given by SE^"^ ~ {dHy/dz)c/ijjK ~ 5jx^Jk'^ — uP' jc^ jion, 
for a wave with the wavevector k along the z-direction. 
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FIG. 3; (color online) Au(lll) film of thickness d embedded 
into an insulator with dielectric constant n. Dark area con- 
tains bulk electrons described by eiiS). Stripes near the film 
surfaces indicate 2DEG of Tamm-Shockley electrons. Arrows 
show directions of E:c for symmetric (s) and anti-symmetric 
(as) modes. 



Thus, the relative magnitude of the two effects is 



(21) 



where in the last identity we made use of Eqs. ([H]) and 
the response functions (fT0| - (fT2l) . To obtain a numerical 
estimate of the ratio in Eq. (|21|) it is sufficient to utilize 
the plasmon spectrum for A = (see Eq. ((29|l with the 
first two terms only). At low frequencies, much less than 
the plasma frequency, w <C 17, we find, nj ^JkP- — {jp jc? sa 
^tanh(^). We thus obtain, 



— ^ q — tanh (nd/2c) . 



(22) 



As seen from this expression, smaller values of spin-orbit 
splitting A favors the role of magnetization. At present 
we do not know the gyromagnetic ratio g of surface elec- 
trons. For g = 2 and fl — 1.4 x lO^^s^^ we estimate 
^ gP£L^ = 0.14 for Au(lll). The effect of the electric 
current is thus considerably stronger than the effect of 
magnetization. For thin films d < c/D, the suppression 

of SEi^^ is even more noticeable. Yet, assuming that in 
some other materials the situation can be different, in 
the following we analyze both contributions. 

For a plasmon wave propagating along the x direction, 
the electric field E(z)e'*'^~*"* obeys the wave equation, 
which inside the film (|z| < d/2) has the form, 



dzJ 



- fc - E = 



(23) 



where the dielectric function e is still given by Eq. ([4|. 
The corresponding wave equation outside the film has 
the same form as Eq. ((23)) with e — > k. Since equa- 
tion (I23p is invariant under spatial inversion transforma- 
tion, z —z, its solutions are either symmetric or anti- 
symmetric under this inversion. The symmetric solution 



can be written as 



Exiz) 



El 



shjKiz) 



\z\ < d/2, 



cosh {Kid/2) ' 
^^^-K2\z\^d/2^ d/2<\z\, 



(24) 



where Ki — ^Jk!^ — euP jc?- and K2 = \Jk'^ — nuP- jc?-. 
The anti-symmetric solution can be written in the same 
manner only with the change, cosh — )■ sinh. 

The magnetic field is found similarly to how it is done 
in Sec. IV and reads 



sEi sinh (Kiz) 
iuj I i^i cosh {K-id/2) ' 



\z\ < d/2, 



sgn(2)^e-'f^-^l^|-'^/^ d/2 < 



(25) 

The boundary conditions relate the discontinuity of the 
electric field to the interface magnetization 



„ „ 4:T:iuj 

E2-El = My, 



(26) 



and the discontinuity of the magnetic field to the surface 
current. 



£Eii&-nh{Kid/2) KE2 _ ^nj^ 



(27) 



cf. Eqs. (fT5|) . As already pointed out in Sec. Ill, the 
right-hand sides of boundary conditions (I26II27P are not 
well defined as the interface current and magnetization 
([9|) have been calculated under the assumption that the 
electric and magnetic fields are uniform across the inter- 
face. Adopting the same approximation as utilized above, 
cf. Eq. ([HI), we write 



Ex{d/2) = liEi+E^), 



(28) 



and similarly for the magnetic field. After straightfor- 
ward transformations we find the dispersion equation for 
symmetric plasmons. 



, tanh|Xi 



4nuj^ £(w)Ktanh 

K1K2 



(29) 



The spectrum of anti-symmetric plasmons is determined 
by the equation which is obtained from Eq. ([29]) by the 
substitution tanh — > coth. Qualitatively, the modifica- 
tion of the spectrum of low-frequency anti-symmetric 
modes is similar to that of symmetric modes though 
quantitatively smaller. Below we discuss symmetric plas- 
mons only. 

The first two terms in Eq. (|29p describe conventional 
plasmons in a metallic film surrounded by a dielectric. 
For d — > cx) and ck ^ to they yield a well-known surface 
plasmon dispersion relation e + n = 0. The third term 
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originates from the electric response of Tamm-Shockley 
surface electrons and contain an intrasubband Drude 
term as well as an intersubband resonant contribution, 
cf . Eq. (fTU]) . Finally, the last term describes the magnetic 
response of surface states, which also features interband 
resonances at = A, see Eq. As can be easily ver- 

ified, the relative magnitude of magnetic and electric 2D 
responses is indeed controlled by the above parameter, 
Eq. jia). 

At low frequencies, a; <JC f2, and the dielectric function 
is e sa juP'. Since wavelengths of interest are small, 
k <^ f^/c, after simple transformations (neglecting the 
last term in Eq. (j29[) . which is typically a good approxi- 
mation), we find a plasmon spectrum in the form. 



KUJ 



cQ tanh ^^ 

2c 



2p^[(u;+z/t)2-A2] 



2 • 



(30) 

Here we assumed that scattering rate l/r is small com- 
pared with frequency u and retained it in the resonant 
intersubband term only. For thin films, d <C c/fl, which 
for Au means that d < lOnm, the electric field is al- 
most uniform across the film, so that oscillations of sur- 
face electron density simply add to the oscillations of the 
charge density in the bulk of the film. 

Fig. 0] illustrates the effect of the interband spin res- 
onance on the plasmon spectrum of thin Au(lll) films. 
Near the point where this spectrum intersects the spin 
resonance frequency, ui = A, there is a significant de- 
crease in the plasmon group velocity, duj/dk, and even 
a narrow region of negative group velocity. In addition, 
at Lu — A the plasmon phase velocity is increased signifi- 
cantly and (in the absence of electron scattering) reaches 
the value of cj^/n,^ characteristic of photon propagation 
in the dielectric. The corresponding increase of the wave- 
length 1/fc due to the resonant interaction with surface 
electrons can be quite significant, as seen from Fig. [4l A 
finite scattering rate 1 /r makes these features less sharp, 
but broadens them, which could in fact help their observa- 
tion. However, the scattering rate should be made rather 
small, which emphasize the use of low temperatures. 



V. PLASMON SPECTRUM OF A METAL FILM 
EMBEDDED IN AN IONIC DIELECTRIC 

Another interesting manifestation of the surface states 
can be achieved if the material surrounding the metal film 
is an ionic crystal. (Here we assume that the presence 
of a dielectric, while possibly changing the value of A, 
does not suppress it significantly.) As well-known, the 
dielectric function of an ionic crystalS^, 



k{uj) 



(31) 



has a pole at the frequency of the transverse optical 
phonon lot and a zero at the frequency of the longitu- 
dinal optical phonon lul. Since the dielectric function is 
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FIG. 4: (color online) The effect of the interband spin res- 
onance on the spectrum of symmetric plasmons in a thin 
Au(lll) film, d = 2 nm, coated with a dielectric with k = 10. 
The inset shows the full frequency range. The main plot shows 
frequencies close to the resonant frequency A — 100 meV. 
The interaction with the surface electrons leads to strong res- 
onant modification of plasmon spectrum close to a; = A. The 
plasmon spectrum is shown for different values of the scat- 
tering rate: l/r = 1 x 10^^ s~^ (blue dashed-dotted line), 
3 X 10^^ s~^ (green dashed line) and 5 x 10^^ s~^ (dotted 
black line). Red dashed and black solid lines (shown also in 
the inset) represent the asymptotes when l/r — > 0. The plas- 
mon group velocity is varied greatly in this frequency range 
and can change sign provided that scattering rate is small 
enough. Note that the plasmon spectrum is bounded by the 
photon line uj — ck/^fn from the left. 



negative within the frequency range, ut < lo < lol^ the 
plasmon (or, more accurately, plasmon-polariton) exci- 
tations normally cannot propagate in this range of fre- 
quencies. This can be verified from Eq. ([23) with A = 0, 
as the first two terms in this equation are negative. 

However, the finite spin-orbit splitting can change 
this situation dramatically if the splitting energy A falls 
within the forbidden frequency band, i.e. lot < A < lul. 
For a Au(lll) film with A « 100 meV such a situa- 
tion can be realized, in particular, for LiH, Hoo = 3.6, 
ivt = 72 meV, ujl = 138 mcV or SiC, Kqo — 6.7, 
LOT = 98 meV, ujl = 117 meV, data from Ref. i23|. Since 
the resonant term in Eq. ([TU)) changes sign below w = A, 
plasmon-polariton modes can propagate in the bandgap 
with a frequency close to the frequency of intersubband 
transitions. 

For thin films d <^ c/Cl the bandgap resonance is at 



= A- 



2dp|, • 



(32) 



The impedance of the coated film can be obtained sim- 
ilarly to the calculations of the preceding Sections. We 
assume here that the thickness of the insulator coating D 
is smaller than the decay length of electromagnetic field 
in the insulator, D <C c/A < 2 fim. The impedance is 
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then found to be 
^H = ^ 1 , J . (33) 

We obtain that the width of the resonance is given by the 
scattering rate and its relative ampHtude is proportional 
to the ratio of the second term in Eq. ([33]) and the scat- 
tering rate. For the metal film of thickness d = 10 nm 
and scattering rate l/r = 10^^ s^^ the relative hight of 
the resonance will be ~ 10^"^. 

Let us emphasize again that the advantage of a set-up 
with a polar dielectric coating lies in the fact that the res- 
onance (|32p is purely due to the intcrsubband transitions 
with no other excitations existing in the same frequency 
domain that could otherwise obscure their observation. 



VI. SUMMARY AND CONCLUSIONS 

Noble metals, such as Au or Bi/Ag alloys, possess 
the unique property that a substantial density of two- 
dimensional surface (Tamm-Shockley) states is present 
when the surface is grown in a particular direction (111) 
of a cubic lattice. This allows the study of the inter- 
play of 2D and bulk 3D electron liquids. Of particular 
significance is the fact that the surface confining poten- 
tial breaks spatial inversion and leads to an intrinsic^ i.e. 
independent of any disorder potential, spin-orbit inter- 
action. The latter results in the formation of spin-split 
surface subbands with the splitting reaching a quite sig- 
nificant magnitude of A ~ 0.1 — 0.2 eV at the Fermi 
level. 

The response functions of the surface states have a 
resonant character for frequencies close to the energy of 
spin-orbit splitting. This resonance allows optical de- 
tection of 2D states, which are otherwise obscured by 
a much large number of bulk electrons. In thin films 
the optical response can be expected to be influenced by 
the presence of surface to a larger extent. Measurements 
of reflection and transmission coefficients in the infrared 
spectrum should be able to reveal the intcrsubband res- 
onances. 

Another phenomenon, which is predicted to bear the 
signature of Tamm-Shockley states, is associated with 
surface plasmons, i.e. collective excitations of electron 
density. Surface plasmons (as well as symmetric plas- 
mons in metal films) are gapless at long wavelengths. 
This ensures that the plasmon energy intersects the en- 
ergy of spin intcrsubband resonances at some wavelength. 
The electric field is then strongly enhanced by the motion 
of 2D electrons which results in significant and detectable 
modifications of the plasmon spectrum. 
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APPENDIX A: SOLUTION OF KINETIC 
EQUATION 

To solve kinetic equation ([7]) written in the frequency 
representation as 

^/p - fp] = -i^p, (Al) 

we first calculate its commutator with 7)n, 

'^[/p, Vn] + Ap(/p - finfpfin) = i[ICp, T]n\, (A2) 

where we utilized that fj^ — 1. Finally, we need one more 
equation, which is found from Eq. (jAl[) 

W77n/p^n + ^Ap[7)n, fp] = ifjnJCpf]^, (A3) 

Formally, the system of three equations (jAmA3[) contains 
three unknowns: /p, [/p,?7n], and finfpfin- Eliminating 
the last two of these unknowns we obtain the solution of 
kinetic equation in the form, Eq. ((8|). 

APPENDIX B: BOUNDARY CONDITIONS 

Let us assume that the amplitude of the electromag- 
netic field acting on surface electrons is small compared 
with the atomic fields (this is a rather weak condition). 
The modification of the electronic wave function can then 
be neglected and the effective electric field in Eq. ([7]) 
acting on surface electrons is written as 

E,^ J dz\iP{z)\''E,{z), (Bl) 

where ipi^) is the wave function describing transverse 
confinement of surface states. Here we adopt the simplest 
approximation that ip{z) is momentum-independent and 
constant across the interface layer of thickness a. Equa- 
tion (|Bip then gives, 

E, = - f dzE^{z). (B2) 
a J a 

The interface layer can now be considered as a metal film 
with the effective dielectric function 

^_ 47re^n2 to* e^A"^ 

m*aLu^ 2uj'^Ppa {to + i/r)^ ~ A'^' 

that is obtained from Eq. (fTU]) : n2 is 2D density of surface 
electrons. Solution of Maxwell's equation for < z < a. 
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which satisfy the condition that E{Q) = Ei and E{a) = 
E2, is 

smh htiZ 

Exiz) = El coshKz + {E2 — El coshKa) — ; , (B4) 

smh na 



used throughout the paper as long as na <^ 1. This 
condition is satisfied as long as 



where n — ^JkP- — euP- As a result, Eq. (jB2p gives the 
effective field, 

„ „ sinh Ka , ^ „ , , cosh Ka — \ , 
Ex = El h {E2 - El coshKa) — • (B5) 



na sinh kg 



max(|w - A|, 1/r) > ^ 10~VeV. (B7) 



This condition reduces to 



Ex — 



El + E2 



(B6) 
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